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We present a good description of recent experimental data on forward di-hadron azimuthal corre- 
lations measured in deuteron-gold collisions at RHIC, where monojet production has been observed. 
Our approach is based on the Color Glass Condensate theory for the small- a; degrees of freedom of 
the nuclear wave function, including the use of non-linear evolution equations with running QCD 
coupling. Our analysis provides further evidence for the presence of saturation effects in RHIC data. 
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Analyses of the angular dependence of two-particle cor- 
relations in hadronic collisions have proven to be an es- 
sential tool for testing the underlying QCD microscopic 
dynamics. Indeed, the presence of well pronounced back- 
to-back azimuthal correlations for high transverse mo- 
mentum particles in proton-proton (p+p) collisions at 
RHIC and Tevatron []| is well understood within the 
standard leading-twist approach to QCD. There, they 
originate from the scattering of two free partons, each 
coUinear to one of the colliding protons, together with 
energy-momentum conservation. In turn, the absence of 
back-to-back correlations observed in gold-gold collisions 
at RHIC 2|, for mid-rapidity pairs, has been interpreted 
as due to the interaction of the outgoing partons with a 
dense and strongly interacting medium produced in such 
collisions, presumably a Quark Gluon Plasma. 

Related measurements have also been performed in 
deuteron-gold (d-|-Au) collisions at RHIC. While for mid- 
rapidity pairs data is consistent with a back-to-back 
structure, recent measurements of di-hadron correlations 
in the deuteron fragmentation region, i.e. at forward ra- 
pidities, display a strong suppression of the away-side 
peak [3|, a phenomenon also referred to as monojet pro- 
duction. The presence of back-to-back correlations at 
mid-rapidity and the fact that total particle multiplici- 
ties decrease at forward rapidities rule out a final state 
interpretation of this novel phenomenon. Rather, they 
suggest an initial state effect or a modification in the pri- 
mary particle production mechanism with respect to the 
standard coUinear factorization approach. 

At high energies, particle production in the forward ra- 
pidity region probe partons carrying a small fraction x of 
the target (nucleus) longitudinal momentum. Moreover, 
the gluon densities within a hadron grow very fast with 
decreasing x. Therefore, and contrary to the assumptions 
of coUinear factorization, gluons in such regime cannot 
be regarded as independent. Rather, they form a highly 
coherent system governed by non-linear phenomena and 
gluon saturation. These qualitative ideas are casted in a 
precise theoretical framework known as the Color Glass 
Condensate (CGC) effective theory of QCD in the high- 
energy limit (see e.g. [J| and references therein). 

There is compelling evidence that gold-gold and d+Au 



collisions at RHIC probe the saturation regime of QCD. 
This claim is supported by the success of saturation mod- 
els p, |6| in the description of the energy, rapidity and 
centrality dependence of the transverse-momentum in- 
tegrated particle multiplicities experimentally measured. 
On the other hand, CGC-based calculations successfully 
predicted first [7,^], and correctly described then [9, li)], 
the measured suppressed production of high transverse- 
momentum forward hadrons in d-|-Au collisions with re- 
spect to p+p collisions [llj, [l^ . Despite the success of 
CGC phenomenology at RHIC outlined above, alterna- 
tive interpretations of data are also possible [13|. Al- 
though some of these approaches invoke dynamical in- 
gredients related to those in the CGC, such as unitarity 
or multiple scatterings, it remains unsettled whether the 
latter is indeed the right approximation to QCD that 
better grasps the relevant dynamics probed at RHIC. 

Disentangling the presently unclear situation requires 
the study of more exclusive observables. Forward di- 
hadron correlations in d+Au collisions are a good test- 
ing ground for models. Indeed, the possibility of monojet 
production as a result of saturation effects was antici- 
pated in IJ, [l5|, albeit at a qualitative level. In this 
letter we shall present a CGC analysis, using the most 
up-to-date theoretical tools available, of the recent data 
on forward di-hadron azimuthal correlations in d+Au col- 
lisions at RHIC by the STAR collaboration, providing a 
good description of the data. We argue that our results 
provide further evidence for the presence of saturation 
effects in RHIC data and also that the CGC effective 
theory approximates well QCD in the saturation regime. 

We are interested in the calculation of the double- 
inclusive hadron production cross section in proton- 
nucleus (p+A) collisions dfjP A^hih2X ^ rpj^^ theoreti- 
cal groundwork needed for this calculation in the CGC 
was presented in [15l4l8l |. We shall denote pi.2i_ and 
2/1,2 = lii(Pi 2/Pr2)/2 the transverse momenta and ra- 
pidities of the produced hadrons, using light-cone co- 
ordinates for four-dimensional vectors p — {p^ ,p^ ,p±). 
We introduce the Feynman-x variables Xi — \pi±\e^^ /y/s, 
with ^/s the collision energy per nucleon. The kinematics 
of di-hadron production is such that the partons in the 
proton (nucleus) wave function that can contribute to the 



cross section carry a fraction of longitudinal momentum 
bounded from below by Xp {xA)^ given by 



Xp — xi + X2 , XA = xi e 



'2yi 



X2 e 



'2^2 



(1) 



The STAR measurements to be analyzed here feature 
^/s — 200 GeV » \Pii_\ and yi ^ 3, yielding Xp^QA and 
XA ~ 10"'^. Thus, the scattering process involves quan- 
tum fluctuations well understood in QCD on the proton 
side, while, on the nuclear side it probes quantum fluctu- 
ations whose non-linear QCD dynamics can be studied. 

First measurements of the coincide nce p robability were 
performed at more central rapidities [1^ |l9[ , less suited 
for exploration of small-x effects. These purely kinematic 
considerations motivated the approach of [15|, which we 
shall follow here. There, the proton is treated as a 
collinear object and described by means of the usual par- 
ton distribution functions, whereas the nucleus target is 
probed in the high parton density regime and is amenable 
to a CGC description. In this picture, the scattering pro- 
cess is initiated by a fast valence quark from the incom- 
ing proton (or deuteron) which splits into a quark and a 
gluon, either before or after scattering off the saturated 
small- a; glue of the target. Importantly, the splitting into 
a quark-gluon pair is described by standard pQCD tech- 
niques and, therefore, features an exact back-to-back cor- 
relation in the transverse plane. The quark-gluon system 
is put on-shell through the interaction with the nucleus, 
as a result of which the quark and gluon also acquire a 
transverse momentum of the order of the saturation scale 
of the nucleus. When that scale, which marks the onset 
of non-linear effects, is comparable to the initial trans- 
verse momenta of the quark and gluon, their intrinsic 
angular correlation is washed out. Finally, the outgoing 
quark and gluon fragment independently into hadrons. 

The forward di-hadron production yield in p+A colli- 
sions reads 



dNP"^^^^^^^ = f dzi f dz2 f dx 

J'-^l J^2 J^ + ^ 



fq/pix,tJ-) 



^^,A^,sX fp^ g^ ^\ D^,^/g{z,,fi)D^,^/^{z2, ^Ji) ,(2) 

where fq/p and D^/i represent the valence quark distri- 
bution of the proton and the fragmentation function of a 
parton i into hadron h respectively. We use the CTEQ6 
NLO quark distributions [l^l and the KKP NLO frag- 
mentation functions ^21| . In the analysis of d-|-Au colli- 
sions we shall neglect all nuclear effects in the deuteron, 
regarded as an incoherent superposition of a proton and 
a neutron. We obtain the neutron valence quark dis- 
tribution from the proton one using isospin symmetry. 
The factorization and fragmentation scales are both cho- 
sen equal to the transverse momentum of the leading 



hadron, which we denote as hadron 1, fi = \pi±\. Fi- 
nally, the cross section for the inclusive production of 
the quark-gluon system in the scattering of a quark with 
momentum xP^ off the nucleus A reads [la : 



flJSfQA^qgX ^^(j^ 



5{xP+-k+-q+) F{£a,^) 



Aq/3 



(3) 



where P, q and k denote the deuteron, quark and gluon 
momenta, respectively. We denoted A = k± + q±, 
z — k^ /xP^ and the longitudinal momentum fraction 
of gluons in the nucleus xa = xi e~'^'^^ /zi +X2 e~'^y^ / Z2- 
Due to parton fragmentation the values of x probed are 
generically higher than Xp and xa defined in Eq. ([1]). For 
the proton, one has 0.4 ^ Xp < x <1. In that region the 
proton wave function is dominated by valence quarks. 
Therefore we have safely neglected the gluon initiated 
processes gA — >• qqX and gA -^ ggX in Eq. ^ . 

Let us review the ingredients in Eq. ([3]). F{x,k±) is 
the unintegrated gluon distribution (ugd) of the target 
nucleus. It is related to the correlator A/" of two light-like 
fundamental Wilson lines in the background of the color 
fields of the nucleus, also referred to as dipole scattering 
amplitude, through a Fourier transform: 



F{x,k^) = 



^ — if*^ I 



(27r)2 



[1-Af{x,r)], (4) 



where r denotes the dipole transverse size. The sec- 
ond line of Eq. (jSj features the /cT-factorization breaking 
term. Its exact evaluation requires the knowledge of four 
and six-point correlators of Wilson lines 15l-ll8|. which 
are difficult to obtain in practice. However, assuming 
that the distribution of color charges in the nucleus is 
Gaussian with a non-local variance, higher-point correla- 
tors can be expressed in terms of the two-point function 
F{x,k±). As shown in 15|, further use of the large-TVc 
limit yields the following expression for /: 

lapi^if^^'^x) = / C?^'7i'0a/3(^>'?i)-P'(^'fc-L-'7±) , (5) 

where ip^a is the amplitude for q —>■ qg splitting and A, 
a and /3 are polarization and helicity indices. We re- 
mark that, even under the Gaussian approximation, the 
cross section ([3]) is still a non-linear function of the un- 
integrated gluon distribution, thereby invalidating kr- 
factorization, namely dN^^^iB^ ^ F (gi dNis^is^ . Fi- 
nally, the delta function in Eq. ^ is due to the use of 
the eikonal approximation. It reflects the fact that in a 
high-energy scattering, the momentum transfer is mainly 
transverse, since the exchange of longitudinal momentum 
is suppressed by powers of the collision energy. 

The CGC is endowed with a set of non-linear evolu- 
tion equations which in the large- A^c limit reduce to the 



Balitsky-Kovchegov (BK) equation 22|, |23|. These equa- 
tions are renormalization group equations for the x evo- 
lution of the unintegrated gluon distribution, and more 
generahy of n-point correlators, in which both linear ra- 
diative processes and non-linear recombination effects. In 
this work, we compute the small-x dynamics of the dipole 
correlator, and hence that of the ugd, by solving the BK 
equation including running coupling corrections (rcBK): 



dU{x, r) 



dVif™"(r,ri,r2)[AA(:E,ri) 



d\n{xo/x) 

+Af{x,r2) - Afix,r) - Af{x,ri)Af{x,r2)] . (6) 



The evolution kernel i^T"^"" 
prescription of 



24| 



is evaluated according to the 
Explicit expressions for the ker- 
nel, together with a detailed discussion on the numerical 
method used to solve the rcBK equation can be found in 
[25j. The only piece of information left to fully complete 
all the ingredients in Eq. ([2]) are the initial conditions 
for the rcBK evolution. This non-perturbative input has 
been constrained by the analysis of single-inclusive for- 
ward hadr on p roduction in d-|-Au collisions at RHIC per- 
formed in [10| using an analogous CGC set up: the start- 
ing point of the rcBK evolution is xq — 0.02, and at that 
value of X the initial saturation scale probed by quarks is 
OgQ = 0.4 GeV^. Here, we simply take over this informa- 
tion. In this respect, the forward di-hadron calculation 
presented here is parameter-free. 

We will now investigate the process dAu — > hih2X, 
with -^3 = 200 GeV. More specifically, we are interested 
in the coincidence probability, an experimental quantity 
measured by both the PHENIX and STAR collaborations 
at RHIC. It is given by CP{A<j)) = Npair{A4>)/Ntrig with 



iVpa^r(A0) 



dN 



dAu^^h\hiX 



Vi,\PiA 



d^pid^P2 



Ni 



trig 



dN 



dAu-^hX 



d^p 



y, p± 



(7) 

and it has the meaning of the probability of, given a trig- 
ger hadron hi in a certain momentum range, produce an 
associated hadron ft,2 in another momentum range and 
with a difference between the azimuthal angles of the two 
particles equal to A0. In order to study the centrality 
dependence of the coincidence probability, we identify 
the centrality averaged initial saturation scale Q%, ex- 
tracted from minimum-bias single-inclusive hadron pro- 
duction data, with the value of Q^q at b = 5.47 fm 26j, 
and use the Woods-Saxon distribution TA{b) to calculate 
the saturation scale at other impact parameters: 



Q%{b) 



TA(5.47fm) 



Qlo = 0.4 GeV^ 



(8) 



Following the experimental analysis by the STAR collab- 
oration we set |pi_l| > 2 GeV, 1 GeV < \p2±\ < \pi±\ and 
2.4 < yi,2 < 4, and require both hadrons to be neutral 
pions. Single-inclusive hadron production, used to nor- 
malize the coincidence probability, is calculated as in [10| . 
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FIG. 1: The coincidence probability as a function of A(/!> 
for central d-|-Au and p-|-p collisions. The preliminary data 
points by the STAR collaboration have been shifted vertically 
by constant amounts. Our results are displayed in solid lines. 

We shall also calculate CP{A(j)) in p-l-p collisions in an 
analogous way, with the initial saturation scale Q^ — 0.2 
GeV^ for the proton. Even though the applicability of 
our approach to the p-l-p case is questionable, the results 
obtained will be useful to interpret the phenomenon of 
monojet production observed in d-|-Au collision. 

Our results for CP(A(/)) in central d-|-Au and in p-fp 
collisions are displayed in Fig. [1] along with preliminary 
data from the STAR collaboration [3[ . Several comments 
are in order. First, the disappearance of the away-side 
peak around A^ ~ tt in d-|-Au collisions exhibited by 
data is quantitatively well described by our CGC calcu- 
lation. In turn our results for p-l-p collisions, where the 
applicability of the CGC formalism is not fully justified 
due to the smallness of the proton saturation scale, agree 
at a qualitative level with the presence of a well defined 
away-side peak in p-l-p collisions. 

Analogous measurements of the coincidence probabil- 
ity in d-|-Au collisions at mid-rapidity, where xa is large, 
also display a clear back-to-back correlation [l2|, |l9| . We 
thus conclude that monojet production is linked to the 
presence of a high gluon density in the target or, equiva- 
lently, to the fact that its saturation scale is comparable 
to the momenta of the produced hadrons. Similar conclu- 
sions have been obtained relying on a saturation model 
[27| (there, although different working assumptions are 
used, the presence of a saturation scale is also the crucial 
ingredient to successfully reproduce data). Our assump- 
tion of independent parton fragmentation in Eq. ([2]) pre- 
vents us from extending our calculation to the near-side 
A(f) ~ region. This would require the use of poorly con- 
strained di-pion fragmentation functions. Moreover, the 
different measurements show little variation in the height 
and width of the near-side peak with varying colliding 
system or centrality, indicating that the near-side peak 
is not sensitive to saturation physics. Finally, since un- 
correlated background has not been extracted from data, 
we have adjusted the overall normalization of data points 
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FIG. 2: The coincidence probability as a function of A(/> for 
different centralities of the d+Au coUisions. 

by subtracting a constant shift. 

To further illustrate how the disappearance of the 
away-side peak is related to the increased gluon density, 
and quantify better the transition between p+p and cen- 
tral d-|-Au collisions, we show in Fig. [2] the centrality 
dependence of CP{A(j)). We predict that the near-side 
peak does not change with centrality, and that the away- 
side peak reappears for less central collisions. This is 
consistent with the fact that peripheral d-f Au collisions 
are essentially p-fp collisions. We also predict that for 
higher hadron transverse momenta the away-side peak 
will reappear, as larger xa values will be probed [15]. 

We conclude that the physics of monojet production is 
due to the interplay between the transverse momenta of 
the produced hadrons and the one acquired during the 
interaction with the nucleus. In the CGC approach pre- 
sented here, the interaction with the nucleus is realized 
in a fully coherent way, and the momentum broadening 
is parametrically controlled by the x-dependent satura- 
tion scale of the nucleus. The latter, in turn, is described 
by means of the rcBK equation. While different inter- 
pretations of data may be possible, presently we are not 
aware of any alternative description of this phenomenon 
that does not invoke saturation effects. On the other 
hand, the simultaneous description of single and double 
inclusive hadron production at RHIC sets a constrain- 
ing test to models. Here we obtain a good description 
of the d-|-Au data without modifying the few parame- 
ters already constrained by the analysis of single- inclusive 
hadron production presented in [1C|. Energy loss effects, 
which were a key ingredient in the alternative explana- 
tions [13] of the suppression of forward hadron yields in 
d-|-Au collisions, are not required in our approach in order 
to obtain a good description of forward di-hadron corre- 
lations. In this respect, the analysis presented here lends 
support to the idea that manifestations of the saturation 
regime of QCD have been observed at RHIC, and that 
such regime is well described by the CGC at its present 
degree of accuracy. 
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